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ABSTRACT 

In  this  paper  the  three  level  formulation  of  a 
variational  multi-scale  (VMS)  large  eddy  simulation 
(LES)  method  for  compressible  flow  computations  is 
extended  for  applications  involving  moving/deforming 
grids.  A  consistent  method  to  improve  the  VMS-LES 
method  by  computing  the  small-scale  Smagorinsky 
constant  dynamically  as  the  flow  develops  (dynamic 
VMS-LES)  is  also  extended  for  dynamic  grid 
applications.  Two  applications  of  VMS-LES  for  the 
simulation  of  separated  flow  over  moving  NACA-0012 
extruded  airfoil  is  then  presented.  The  first  application 
involves  a  qualitative  simulation  exploring  the  Knoller- 
Betz  effect  of  the  heaving  airfoil  at  high  Strouhal  number. 
The  second  application  is  that  of  the  pitching  airfoil 
undergoing  dynamic  stall.  The  results  predicted  by  the 
dynamic  VMS-LES  method  are  compared  to  those 
obtained  with  other  turbulence  models  and  to 
experimental  data  and  it  is  found  that  the  dynamic  VMS- 
LES  performs  better  than  the  other  considered  static  and 
dynamic  LES  models. 


1.  INTRODUCTION 

Since  its  introduction  the  VMS-LES  method  (Hughes 
et.al,  2000,  2001)  has  been  successfully  applied  to 
simulate  many  practically  relevant  turbulent  flow 
computations.  In  this  approach  the  Navier-Stokes  (NS) 
equations  are  treated  by  a  variational  projection  instead 
of  the  classical  filtering  methodology.  As  clarified  later 
in  Collis,  2001,  in  VMS-LES,  spatial  scales  are  separated 
a  priori  and  systems  of  equations,  each  representing  the 
dynamics  of  one  of  the  separated  scales,  are  derived. 
Thus,  this  formalism  offers  the  freedom  to  choose 
different  closure  models  for  equations  representing 
different  scales.  The  system  of  equations  can  then  be 
recombined  to  obtain  the  final  variational  form  of  the  NS 
equations  appended  with  closure  models.  Typically  in 
VMS-LES,  the  effect  of  the  unresolved  scales  is  modeled 
only  in  the  equations  representing  the  smallest  resolved 
scales.  Therefore,  no  energy  is  directly  extracted  from 
the  large  structures  in  the  flow.  These  class  of  methods 
have  shown  great  promise  in  simulation  of  different 
turbulent  flows,  such  as  wall  bounded  flows  (Hughes 
et.al,  2001)  ,  separated  turbulent  flows  (Koobus  and 
Farhat,  2004,  Rajasekharan  et.al.  2007,  2008)  ,  turbulent 
flow  control  (Ramakrishnan  and  Collis,  2002)  and 


simulation  of  bypass  transition  (Calo,  2004)  to  name  a 
few. 

Many  variants  of  the  VMS-LES  method  have  been 
developed  over  the  last  few  years.  Initial  attention 
focused  on  incompressible  flows,  regular  grids  and 
spectral  discretizations  wherein  the  a  priori  separation  of 
scales  was  achieved  via  a  frequency  cutoff.  For  finite 
element  approximations,  a  scale  separation  algorithm 
based  on  hierarchical  bases  was  presented  in  Jansen  and 
Martinez,  2002.  For  compressible  turbulent  flows,  a 
more  general  approach  for  discretizing  the  two-level 
VMS-LES  formulation  on  unstructured  finite  volume 
(FV)  as  well  as  finite  element  (FE)  meshes  was  proposed 
in  Koobus  and  Farhat,  2004.  Later  in  three-level 
formalism  in  tandem  to  the  approach  in  Farhat  et.al., 
2006  for  compressible  flows  was  derived  and  was  shown 
to  introduce  requirement  of  closure  terms  in  the 
continuity  equation  as-well. 

VMS-LES  involves  two  stages,  first  is  a  formalism  stage 
as  described  in  the  first  two  paragraphs  and  the  second  is 
a  closure  stage.  In  most  previous  works  on  the  VMS- 
LES  method,  the  closure  was  achieved  by  extracting 
energy  from  the  small  resolved  scales  by  the 
Smagorinsky  eddy  viscosity  model  (Smagorinsky,  1961) 
with  a  constant  coefficient.  In  Farhat  et.al,  2006,  a 
dynamic  VMS-LES  method  was  constructed  by 
computing  this  coefficient  dynamically  (Holmen  et.al, 
2004)  using  the  same  procedure  as  in  the  original 
dynamic  subgrid-scale  eddy  viscosity  model  (Germano 
et.al.,  1991).  In  that  study,  it  was  concluded  that  the 
dynamic  multi-scale  model  was  less  sensitive  to  the  scale 
partition  when  compared  to  its  static  counterpart.  In 
Gravemier  et.al,  2004,  an  iterative  procedure  based  on 
the  concept  of  residual-free  bubbles  (Franca  et.al,  1998) 
was  proposed  for  adjusting  dynamically  the  value  of  the 
coefficient  of  the  Smagorinsky  eddy  viscosity  model.  An 
alternative  approach  for  constructing  a  dynamic  VMS- 
LES  method  that  is  based  on  a  simple  variational 
interpretation  of  Germano’s  identity  was  proposed  in 
Farhat  et.al,  2006.  The  salient  features  of  this  alternative 
approach  are  numerous.  They  include:  (a)  consistency 
with  the  variational  aspect  of  the  VMS-LES  method  as 
opposed  to  the  straightforward  procedure  adopted  in  for 
incorporating  a  dynamic  subgrid-scale  eddy  viscosity 
model,  (b)  consistency  of  the  value  of  the  dynamically 
computed  Smagorinsky  constant  with  the  adopted 
discretization  method  (similar  to  the  vector-level  identity 
presented  in  Oberai  and  Wanderer  2005  and  Morinishi 
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and  Vasilyev  2002),  and  (c)  incorporation  of  the  time- 
history  and  diffusive  effects  by  incorporation  of  the 
inertial  and  viscous  terms  in  the  computation  of  the 
dynamic  Smagorinsky  constant. 

With  the  increase  in  computational  power  over  the  years, 
and  the  development  of  near  wall  models,  it  is  now 
feasible  to  apply  LES  simulations  to  study  flow  over 
small  sized  fixed  and  flapping  wings  in  the  low  to 
moderate  speed  aerodynamics  regime.  This  development 
has  opened  up  an  opportunity  to  successfully  study  flow 
over  Micro  Air  Vehicles  (MAV).  It  is  thus  worthwhile  in 
exploring  the  efficacy  of  the  VMS-LES  method  in 
capturing  the  unsteady  aerodynamic  flows  over  moving 
boundaries.  This  paper  is  hence  an  attempt  to  study  the 
applicability  of  the  VMS-LES  method  for  flows  on 
moving  and  deforming  grids.  In  this  respect,  this  paper  is 
essentially  a  sequel  to  the  previous  works  done  in  . 

This  paper  is  organized  as  follows.  Section  summarizes 
the  arbitrary  Lagrangian-Eulerian  (ALE)  governing 
equations  and  the  FV/FE  formulation.  Section  describes 
the  VMS-LES  formulation  of  compressible  turbulent 
flows  developed  in  Koobus  and  Farhat,  2004  and  re¬ 
derived  in  Farhat  et.al.  2006.  In  particular,  its 
agglomeration-based  procedure  for  separating  a  priori 
the  scales  on  unstructured  meshes  is  highlighted.  In 
Section,  an  overview  of  a  dynamic  extension  of  the 
VMS-LES  based  on  the  variational  interpretation  of 
Germano's  identity  is  presented.  Next  in  Section,  the 
specifics  of  the  spatial  and  temporal  discretizations 
adopted  in  this  work  are  discussed.  In  Section,  the  results 
from  the  application  of  the  dynamic  VMS-LES  method 
to  the  simulation  a  heaving  and  pitching  extruded 
NACA-0012  airfoil  is  presented  and  the  results  are 
discussed.  Finally,  Section  concludes  this  paper  with  a 
summary  of  the  findings  and  reflections  on  future 
improvements  and  challenges. 

2.  ARBITRARY  L AGRAN GIAN  EULERIAN  (ALE) 
FORMULATION  OF  THE  GOVERNING 
EQUATIONS  AND  ITS  SEMI-DISCRETIZATION 

Then  the  ALE  strong  conservation  form  of  set  of  mass, 
momentum  and  energy  equations  governing  a 
compressible  fluid  can  be  written  down  as 

M  +  Jv.(p(  u-w))  =0 

'  d(^U)  +  JV.(/>u®(u-w))  =  -VP  +  V.< t  (1) 

^4^-  +  JV.[(E(u  -  w)  +  Pu)]  =  V.(<ru)  +  V.(AVT) 


Without  any  loss  of  generality  as  far  as  issues  pertaining 
to  unstructured  meshes,  the  following  spatial 
discretization  is  assumed: 

•  The  bounded  flow  domain  is  discretized  by  a 
tetrahedral  mesh  from  which  a  dual  mesh  defined  by 
cells  or  control  volumes  is  derived. 

•  The  convective  fluxes  are  approximated  by  a  finite 
volume  (FV)  method 

•  The  diffusive  fluxes  are  approximated  by  a  piecewise 
linear  finite  element  (FE)  method 

Now,  integrating  Eq.  (1)  on  an  elementary  volume  of  the 
material/reference  space  and  transforming  back  to  the 
physical  space,  the  following  mixed  FE/FV  formulation 
is  obtained 

IA(Xi,  W,w)  =0 

B(^,^,W,w)  =0  (2) 

C{Xi,$u  W,w)  =0 

where 

A(A;,W,w)  =  [  ^XidQ  +  /  p(u-w).n,*;dr 

Jn(t)  dt  Jdsiiit)  v 

)=/  ^XidQ  +  f  p  u  (g>  (u  —  w)  nXi  dT 

Jn(t)  dt  JdUi(t) 

+  [  PnXidT  +  [  crV^idn  (3) 

Jdsii(t)  Jn(t) 

=  /  ^rXidU.  +  [  (E(u-Mv)  +  Pu).nXidT 
Jn(t)  at  JdQi(t) 

+  /  +  [  A  VT.V<Mft 

.  Jn(t)  Jn(t) 

and 

W  =  (p,u,T)* 


3.  VMS-LES  FORMULATION  OF  THE  ALE 
GOVERNING  EQUATIONS 

The  decomposition  a  priori  of  these  spaces  into  large 
resolved  scales,  small  resolved  scales,  and  unresolved 
scales  (a.l.a.  Collis)  can  be  written  as 


VFv 

-  Vfy  ©  V'fv  ®  Vfv 

(4) 

Vfe 

=  Vfe  ®  V'fe  ®  Vfe 

(5) 

where  the  notation  64  ’  ”,  and  “A”  designates  the  large 
resolved,  small  resolved,  and  unresolved  scales, 
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Plugging  these  decompositions  in  Eq.  (2)  a  coupled  set 
of  equations  for  large,  small  and  unresolved  scales  are 


obtained.  Neglecting  the  effect  of  unresolved  scales  on 
large  scales  and  recombining  the  large  and  small  scale 
equations  the  following  set  of  discretized  (subscript  ’h’) 
equations  is  achieved 


'  A(Xih,  Wh,w)  +  A'(X',  W',W,w)  +  RC(X',  W,w) 

+  Cc(X(,  W,W,w)  =0 

B^^W/^w)  +  B'(A",$',  W',  W.w)  +  Rm(A'/,  W,  w) 

+  Cm(^/,$',W,W,w)  =0 

C(Xih,$ih,  Wh)  +  C\X',%W\W,w)  +  Re(X;,*’,W,w) 

+  C.(AJ,*5,  W,W,w)  =0 

where  it  is  noted  that  unlike  the  Favre-averaged 

formulation  of  the  governing  equations,  the  continuity 
equations  requires  modeling.  Choosing  to  ignore  this 
modeling,  and  applying  Smagorinsky  type  closures 

yields 

IA(Xih,  W*,w)  =0 

B(A^,*fc,Wft,w)  +  =0  (7) 

with, 

ms(4>;,w',(c;a;)2)  =  ^t»(c'a;)2|s'|(2sa  -  |s'tJfc<sy)v^'dn 

(8) 

A/„($',W',(C;A;)2,Prt)  =  /  Cp^C‘DA')2lS'l  vr.v-j>;rfn 
'  '  J£l  Pv  t 


3.1  ^  PRIORI  SEPARATION  OF  SCALES  ON 
UNSTRUCTURED  ENTITIES 

A  priori  separation  of  the  scales  according  to 
VpE/FVh  =  VpE/FVh  ®  V*FEfFVh 

is  achieved  here  by  defining  a  projection 

wh=W 

For  this,  control  volumes  and  their  neighbors  in  the  dual 
mesh  are  agglomerated  to  form  macro-cells.  Then,  as 
proposed  in  Koobus  and  Farhat,  2004,  the  large-scale 
component  is 


wfc  =  £**w*  =  E** 

k  k 

with, 


(Y.v°‘(Cj)wA 

jeh 


E  Vo1(cj) 

\  je/fc 


=  E  ®kHWk, 


(9) 


A  priori  separation  of  scales  also  raises  the  question  of 
which  boundary  conditions  to  prescribe  for  the  large  and 
small  resolved  scales.  This  issue  is  addressed  here  by 
requiring  that  the  large  resolved  scales  satisfy  the  no-slip 
(Dirichlet)  boundary  conditions  and  by  allowing  the 
small  scales  to  vanish  at  the  wall.  This  is  consistent  with 
the  observation  that  the  small  eddies  (fluctuations) 
disappear  at  the  wall.  Hence,  the  following  boundary 
restriction  is  incorporated  into  the  definition  of  the  large 
and  small  resolved  scales 


I^FVh/FEh  -  {U/i  •  Vh  -  'E$khVkh,Uh |rwaII  -  u|ru,oH} 
V FVh/FEh  =  { vh  =  _  ®kh)Vkh,  U^r^ii  =  0} 


(11) 


The  above  no-slip  condition  is  enforced  naturally  by 
choosing  not  to  agglomerate  the  first  layer  of  control 
volumes  adjacent  to  the  wall.  This  approach  avoids 
discontinuities  of  the  large  scale  velocities  at  the  wall  in 
comparison  to  strongly  prescribing  the  no-slip  conditions. 
In  this  way  the  Gibbs  phenomenon  is  averted. 

It  should  also  be  mentioned  that  for  computations 
involving  wall-functions  (as  in  moderate  to  high 
Reynolds  number  computations  presented  in  this  paper), 
the  VMS-LES  approach  does  not  need  any  extra 
treatments  at  the  displaced  wall  boundary.  This  is 
because,  the  wall-function  introduces  a  Dirichlet-to- 
Neumann  transformation  which  results  in  imposing  a 
Neumann  boundary  condition  at  the  displaced  boundary. 

4.  DYNAMIC  EXTENSIONS  OF  THE 
VARIATIONAL  MULTISCALE  METHOD 

In  the  previous  section  a  static  VMS-LES  formulation 
for  simulation  of  compressible  turbulent  flows  was 
presented.  One  of  short  comes  of  this  approach  is  that  the 
Smagorinsky  constant  is  fixed.  Hence,  as  it  is,  this  static 
form  of  the  VMS-LES  model  is  not  adaptable  to  the 
different  flow  features  in  the  domain.  Also  this  method 
generates  turbulent  viscosity  even  in  the  laminar  regions 
of  the  flow  domain.  To  overcome  these  shortfalls, 
dynamic  versions  of  the  VMS-LES  method  have  been 
put  forth.  This  section  briefly  describes  the  residual 
based  dynamic  procedure  in  Farhat  et.al,  2006. 

4.1  A  RESIDUAL  BASED  DYNAMIC  VMS-LES 
MODEL  USING  THE  VARIATIONAL 
ANALOGUE  OF  GERMANO’S  IDENTITY 


T* 

£  Voi(Cj)^  Jh 
j€h 


Consider  two  nested  meshes  with  element  sizes  hi  and 
(1°)  h2>hi  respectively,  and  corresponding  finite  dimensional 

subspaces 

C  VFEhl 


3.2  WALL  BOUNDARY  TREATMENS  OF  LARGE 
AND  SMALL  SCALE  VARIABLES 


3 


Substituting  h  by  hx  in  the  second  and  third  of  Eqs.  (6) 
the  momentum  and  continuity  equations  for  each  node  i 
is  obtained.  Similarly,  discretizing  Eqs. (6)  on  the  coarser 
mesh,  but  using  the  prolongated  and  their  associated 
nodal  values  on  the  coarse  mesh,  the  corresponding 
momentum  and  continuity  equations  for  each  node  i  is 
obtained.  Subtracting  the  first  of  this  equations  from  the 
second  and  defining 

-  B(4^vW,,w) 

=  Ms($'ki,w;i,C'Af)  - 

=  MH(Vihi,W'hl,C'A?,Prt)  - 

(14) 

The  above  two  equations  relate  the  resolved  turbulent 
scales  to  the  subgrid-scales  at  two  different  levels.  Hence, 
they  can  be  interpreted  as  a  variational  form  of 
Germano's  identity,  from  which  the  parameters 
Smagorinsky  constant  and  the  turbulent  Prandtl  number 
can  be  dynamically  computed.  Choosing  subspaces 

VFV/FEh2  =  VFV/FEh2  ©  yFV/FEh2  =  VpV/FE ®  [v]  FV/FEh,  (15) 

are  constructed  using  a  second-level  agglomeration 
generated  by  applying  recursive  agglomerations. 

5.  NUMERICAL  SOLUTION 

The  static  and  dynamic  versions  of  the  VMS-LES 
method  presented  in  the  previous  sections  were 
implemented  in  the  AERO-F  module  of  the  nonlinear 
aeroelastic  simulation  platform  AERO-F  (Farhat  et.al, 
2003).  AERO-F  is  a  domain-decomposition-based, 
massively  parallel,  unstructured,  three-dimensional, 
Navier-Stokes  compressible  flow  solver  with  built  in 
moving  mesh  capabilities.  It  blends  an  upwind  scheme 
for  the  convective  fluxes  based  on  Roe’s  approximate 
Riemann  solver  with  a  PI  finite  element  Galerkin 
approximation  of  the  diffusive  fluxes  and  source  terms. 
Higher-order  spatial  accuracy  is  achieved  through  the  use 
of  a  multidimensional  piecewise  linear  reconstruction 
that  follows  the  principle  of  the  Monotonic  Upwind 
Scheme  for  Conservative  Laws  (MUSCL)  (van  Leer, 
1979).  Stabilization  is  obtained  by  a  numerical  diffusion 
that  is  based  on  sixth-order  derivatives  (Debiez  and 
Dervieux,  1999.  The  associated  numerical  dissipation 
has  a  localized  effect  on  high  frequencies,  which  reduces 
its  interaction  with  subgrid-scale  modeling  in  turbulent 
flow  simulations.  Within  this  fluid  module,  RANS,  LES 
and  Hybrid  turbulence  models  can  be  coupled  with 
Reichardfs  wall  function  (Hinze,  1959).  AERO-F  has  the 
capability  of  performing  time-integration  with  a  three- 
point  backward-difference  implicit  scheme  as  well  as  a 
fourth-order  Runge-Kutta  explicit  scheme.  In  the  low 
speed  limit,  it  resorts  to  Roe-Turkel  preconditioning 


(Viozat,  1997;  Turkel,  1987)  of  dissipation  terms  to 
overcome  the  numerical  difficulties  that  are  usually 
encountered  by  compressible  flow  solvers. 

6.  VMS-LES  OF  TURBULENT  FLOWS  ON 
MOVING  GRIDS 

The  static  and  dynamic  versions  of  the  VMS-LES 
method  for  moving  grids  presented  in  the  previous 
sections  are  applied  to  heaving  and  pitching  of  a  NACA- 
0012  extruded  airfoil.  These  flow  configurations  are 
useful  in  the  analysis  and  design  of  MAV  wing  sections. 
All  simulations  presented  in  this  section  utilize  an 
arbitrary  Lagrangian/Eulerian  form  of  Reichardfs  wall 
law  (Tran  et.al,  1998)  and  an  implicit  second-order  time¬ 
stepping  scheme  that  is  provable  second-order  time- 
accurate  on  moving  meshes  and  obeys  its  discrete 
geometric  conservation  law  (DGCL)  (Geuzaine  et.al, 
2003).  The  results  obtained  using  the  VMS-LES  methods 
are  also  compared  with  the  classical  constant  coefficient 
Smagorinsky  LES  model  (Smag  LES),  the  classical 
Favre  averaged  Germano  dynamic  LES  model  (dynamic 
LES,  Camarri  et.al.  2002)  and  to  experimental  data 
wherever  possible. 

6.1  VMS-LES  SIMULATION  OF  THE  KNOLLER- 
BETZ  EFFECT 

The  Knoller-Betz  effect  symbolizes  the  ability  of 
sinusoidally  plunging  airfoil  to  produce  thrust.  This  is  an 
important  consideration  for  the  design  of  propulsion 
systems  of  MAVs  due  to  the  excessive  weight  associated 
with  traditional  battery  driven  fixed  wing  propulsive 
systems.  Various  experimental  and  numerical 
investigations  of  this  phenomenon  have  been  carried  out 
over  the  years  to  understand  the  evolution  of  thrust 
indicative  wake  structures  .  It  was  observed  both 
numerically  and  experimentally  in  Jones  et.al,  1998  that 
the  wake  patterns  at  Strouhal  numbers  greater  than  1.0 
were  non-symmetric,  indicating  that  viscous  effects/flow 
separation  may  be  important  in  understanding  the  wake 
evolution.  This  sub-section  delves  into  the  study 
qualitatively  the  evolution  of  these  wake  structures  using 
the  static  VMS-LES  method. 

To  this  end,  a  qualitative  simulation  of  a  plunging 
extruded  NACA-0012  airfoil  at  a  Reynolds  number  of 
61,000  based  on  the  chord  length  of  1.22m  is  performed 
here.  The  airfoil  is  plunged  with  a  reduced  frequency  of 
k=12.3  and  a  non-dimensional  plunge  amplitude  of 
h=0.12  (resulting  in  a  Strouhal  number  of  1.476).  The 
computational  domain  is  delimited  by  a  circular  cylinder 
of  diameter  20  times  chord  length.  Symmetric  boundary 
conditions  are  imposed  at  the  two  wing  tips.  The  mesh 
has  14,450,920  tetrahedra  and  2,516,690  nodes. 
Preferential  mesh  refinement  is  introduced  in  the  regions 
of  the  separated  wake  as  shown  in  Fig.  1.  below. 
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Fig  1:  Refinement  around  the  NACA-0012  airfoil 
As  seen  in  Fig.  2  the  final  wake  pattern  shows  the  dual 
mode  vortex  street  as  observed  in  ,  which  is  indicative  of 
thrust  and  lift  as  seen  in  the  plot  of  averaged  lift  and  drag 
(Fig.  3). 


Fig  2:  Contour  plots  and  wake  pattern  showing  the  dual 
mode  vortex  street  of  the  plunging  NACA-0012  airfoil 
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Fig  3:  Time  averaged  lift  and  drag  of  plunging  NACA- 
0012  airfoil  using  VMS-LES 


6.2  VMS-LES  simulation  of  the  dynamic  stall  of  a 
NACA-0012  extruded  airfoil 

Many  low  speed  aerodynamics  applications  involve 
highly  unsteady  hovering  flight  regimes  that  call  for  the 
delay  of  aerodynamic  stall.  Dynamic  stall  is  an  unsteady 
phenomenon  where  pitching  a  wing  beyond  its  static  stall 
angle  results  in  maximum  values  of  lift,  drag  and 
moment  coefficients  that  far  exceed  those  achieved  in  the 
static  case.  This  phenomenon  is  characterized  by  a 
massive  flow  separation  and  a  reattachment  process  that 
results  in  a  highly  nonlinear  hysteresis. 

In  order  to  test  the  performance  of  the  different  LES 
models  at  predicting  large  scale  unsteady  effects  and  to 
gain  insight  into  the  phenomenon  of  dynamic  stall,  the 
flow  over  a  pitching  NACA-0012  wing  in  deep  stall 
regime  is  pursued  here.  The  obtained  numerical  results 
are  compared  to  the  experimental  data  reported  in 
McCroskey,  1976.  The  considered  unswept  wing  section 
has  a  span  of  lm  and  a  chord  length  of  1.22m.  It  is 
sinusoidally  pitched  about  the  quarter  chord  so  that  the 
angle  of  attack  is  varied  as 

a  =  15°  +  10°sm(utf) 
where  the  reduced  frequency  is 

k  =  uj  x  Cr/Uoo  =  0.15 

The  free-stream  Mach  number  is  set  to  0.09  and  the 
Reynolds  number  based  on  the  chord  length  is  set  to 
Re=2,600,000. 


Fig  4:  Computational  mesh  for  the  pitching  NACA-0012 
wing  section 

The  computational  domain  is  delimited  by  a  circular 
cylinder  of  diameter  20  times  the  chord  length. 
Symmetric  boundary  conditions  are  imposed  at  the  two 
wing  tips.  The  mesh  has  8,557,236  tetrahedra  and 
1,504,611  nodes.  Preferential  mesh  refinement  is 
introduced  in  the  regions  of  the  separated  wake  as  shown 
in  Fig.4.  For  this  flow,  the  estimated  Taylor  microscale  is 
2.98e-3m,  the  estimated  Kolmogorov  microscale  is 
1.94e-4m  and  the  estimated  Kolmogorov  time  scale  is 
2.6e-5s.  The  characteristic  length  of  the  smallest  grid 
spacing  near  the  wing  and  the  largest  grid  spacing  in  the 
far  field  are  2e-4  m  and  9.2e-2m,  respectively.  The  time- 
step  for  the  simulation  was  set  to  2e-4s. 
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The  number  of  layers  is  set  to  one  in  the  first  and  second- 
level  agglomerations  and  the  depth  of  the  second-level 
agglomeration  is  set  to  one  for  the  dynamic  VMS-LES. 
The  lift  and  pitching  moment  hysteresis  loops  for  the 
deep  stall  considered  here  is  then  computed  using  the 
four  different  LES  models,  namely  the  Smag  LES, 
dynamic  LES,  VMS-LES  and  dynamic  VMS-LES. 
These  loops  are  plotted  in  Fig. 5,  along  with  the 
experimental  data  reported  in  McCroskey,  1976.  It  is 
observed  that  while  all  LES  models  deliver  good  lift 
hysteresis  predictions,  the  dynamic  VMS-LES  model 
delivers  the  best  predictions  for  the  pitching  moment 
coefficient.  Furthermore,  it  is  noted  that  the  dynamic 
form  of  VMS-LES  significantly  improves  the 
performance  of  its  static  counterpart  in  the  prediction  of 
hysterisis. 


l-'l  vs  * 


Fig  5:  Refinement  around  the  NACA-0012  airfoil 


7  CONCLUSIONS 

The  static  and  dynamic  VMS-LES  method  developed  in 
has  been  applied  to  simulate  the  Knoller-Betz  effect  of  a 
heaving  airfoil  and  to  simulate  the  dynamic  stall  of  a 
pitching  airfoil.  The  obtained  numerical  results  suggest 
the  superior  performance  of  VMS-LES  as  a  viable  LES 
model  for  applications  that  involve  turbulent  flow 
computations  over  moving/deforming  objects.  The 
obtained  numerical  results  for  the  dynamic  stall  of  the 
pitching  airfoil  suggest  that  this  dynamic  VMS-LES 
method  tends  to  deliver  more  accurate  results  than  the 
static  VMS-LES  and  the  classical  LES  and  dynamic  LES 
methods.  In  particular,  it  is  found  that  the  dynamic  VMS- 
LES  method  significantly  improves  the  performance  of 
its  static  counterpart  for  the  prediction  of 
lift/ drag/moment  characteristics. 


Predicting  accurate  space-time  correlation  is  an 
important  feature  that  LES  turbulence  models  need  to 
satisfy.  The  dynamic  version  of  the  VMS-LES  method 
developed  here  is  an  attempt  to  bring  in  the  effect  of  time 
correlations  into  LES  simulation,  and  hence  it  would  be 
important  to  study  if  this  model  gives  good  space-time 
correlation.  It  would  also  be  of  interest  to  study  turbulent 
flow  simulations  using  the  VMS-LES  method  to 
characterize  the  flow  over  flexible  structures,  there  by 
working  in  the  full  realm  of  fluid- structure  interaction. 


REFERENCES 

Calo  VM.,  2004:  Residual-based  multi-scale  turbulence 
modeling,  finite  volume  simulations  of  bypass 
transition,  PhD.  Dissertation,  Stanford  University. 

Camarri  S,  Salvetti  MV,  Koobus  B,  Dervieux  A.,  2002: 
Large-eddy  simulation  of  a  bluff  body  flow  on 
unstructured  grids,  Int.  J.  Numer.  Meths.  Fluids 
40:1431-1460. 

Collis  SS.,  2001:  Monitoring  unresolved  scales  in  multi¬ 
scale  turbulence  modeling,  Phys.  Fluids;  1800-1806. 

Debiez  C,  Dervieux  A.,  1999:  Mixed  element  volume 
MUSCL  methods  with  weak  viscosity  for  steady  and 
unsteady  flow  calculation,  Comput.  &  Fluids,  29:89- 
118 

Dervieux  A.,  1985:  Steady  Euler  simulations  using 
unstructured  meshes,  Von  Karman  Institute  Lecture 
Series. 

Dubois  T,  Jauberteau  F,  Temam  R.,  1999:  Dynamic 
multilevel  methods  and  the  numerical  simulation  of 
turbulence,  Cambridge  University  Press. 

Farhat  C,  Rajasekharan  A,  Koobus  B.,  2006:  A  dynamic 
variational  multi-scale  method  for  large  eddy 
simulations  on  unstructured  meshes,  Comput.  Meths. 
Appl.  Mech.  Engrg;  195:1667-1691 

Farhat  C,  Geuzaine  P,  Brown  G.,  2003:  Application  of  a 
three-field  nonlinear  fluid-structure  formulation  to 
the  prediction  of  the  aeroelastic  parameters  of  an  F- 
16  fighter,  Comput.  &  Fluids;  32:3-29. 

Franca  L,  Farhat  C,  Lesoinne  M,  Russo  A.,  1998: 
Unusual  stabilized  finite  element  methods  and 
residual-free  bubbles,  Intemat.  J.  Numer.  Meths. 
Fluids;  27:159-168. 

Germano  M,  Piomelli  U,  Moin  P,  Cabot  WH.,  1991:  A 
dynamic  subgrid-scale  eddy  viscosity  model,  Phys. 
Fluids,  A3: 1760- 1765. 

Geuzaine  P,  Brown  G,  Harris  C,  Farhat  C.,  2003: 
Aeroelastic  dynamic  analysis  of  a  full  F-16 
configuration  for  various  flight  conditions,  AIAA  J. 
41:363-371. 

Geuzaine  P,  Grandmont  C,  Farhat  C.,  2003:  Design  and 
analysis  of  ALE  schemes  with  provable  second- 
order  accuracy  for  inviscid  and  viscous  flow 
simulations,  J.  Comp.  Phys.  191:206-227. 


6 


Ghosal  S.,  Lund  TS,  Moin  P,  Akselvoll  K.  1995:  A 
dynamic  localization  model  for  large-eddy 
simulation  of  turbulent  flows,  J.  Fluid  Mech., 
286:229-255. 

Gravemeier  V,  Wall  WA,  Ramm  E.,  2004:  A  three  level 
finite  element  method  for  the  instationary 
incompressible  Navier-Stokes  equations,  Comput. 
Meths.  Appl.  Mech.  Engrg.  193:1323-1366. 

Gravemier  V.,2006:  Scale- separating  operators  for 
variational  multi-scale  large  eddy  simulation  of 
turbulent  flows,  J.  Comp.  Phys,  212:400-435. 

Hinze  JO.,  1959:  Turbulence,  McGraw-Hill,  NewYork. 

Holmen  J,  Hughes  TJR,  Oberai  AA,  Wells  GN.,  2004: 
Sensitivity  of  the  scale  partition  for  variational 
multi-scale  LES  of  channel  flow,  Phys.  Fluids, 
16(3):824-827. 

Hughes  TJR,  Mazzei  L,  Jansen  KE.,  2000:  Large  eddy 
simulation  and  the  variational  multi-scale  method, 
Comput.  Vis.  Sci,  3-47. 

Hughes  TJR,  Mazzei  L,  Oberai  AA,  Wray  AA.,  2001: 
The  multi-scale  formulation  of  large  eddy 
simulation:  decay  of  homogeneous  isotropic 
turbulence,  Phys.  Fluids,  13:505-512. 

Hughes  TJR,  Oberai  AA  ,  Mazzei,  2001:  L.  Large  eddy 
simulation  of  turbulent  channel  flows  by  the 
variational  multi-scale  method,  Phys.  Fluids 
13:1784-1799. 

Jansen  KE,  Tejada-Martmez  AE.  2002:  An  evaluation  of 
the  hierarchical  basis  in  variational  multi-scale  LES, 
AIAA  Paper  No.  2002-0283. 

Jones  KD,  Dohring  CM,  Platzer  MF,  1998:  Experimental 
and  computational  investigation  of  the  Knoller-Betz 
effect  AIAA  J.;  36(7):  1240-1246. 

Koobus  B,  Farhat  C.,  2004:  A  variational  multi-scale 
method  for  the  large  eddy  simulation  of 
compressible  turbulent  flows  on  unstructured 
meshes  -  application  to  vortex  shedding,  Comput. 
Meths.  Appl.  Mech.  Engrg.;  193:1367-1384. 

Lallemand  MH,  Steve  H,  Dervieux  A.,  1992: 

Unstructured  multigridding  by  volume 
agglomeration:  current  status,  Comput.  &  Fluids; 
21:397-433. 

Mallat  S.,  1999:  A  wavelet  tour  of  signal  processing, 
Academic  Press. 

McCroskey  WJ,  Carr  LW,  McAlister  KW,  1976: 
Dynamic  stall  experiments  on  oscillating  airfoils, 
AIAA  J.;  14:57-63. 

Morinishi  Y,  Vasilyev  OV.,  2002:  Vector  level  identity 
for  dynamic  subgrid  scale  modeling  in  large  eddy 
simulation,  Phys.  Fluids;  14:3616-3623. 

Oberai  AA,  Wanderer  J.,  2005:  A  dynamic  approach  for 
evaluating  parameters  in  a  numerical  method,  Intemt. 
J.  Numer.  Meths.  Engrg.;  62:50-71. 

Oberai  AA,  Wanderer  J.,  2005:  Variational  formulation 
of  the  Germano  identity  for  the  navier-stokes 
equations,  J.  Turb.;  6:1-17. 


Rajasekharan  A,  Farhat  C,  and  Bou-Mosleh  C.,  2007: 
Application  of  a  dynamic  variational  multi-scale 
method  to  the  LES  of  separated  turbulent  flows, 
AIAA  Paper  2007-726. 

Rajasekharan  A.,  2008:  Variationally  consistent  multi¬ 
scale  formulations  and  ALE  time  integrators  for 
Large  Eddy  Simulation  of  turbulent  flows  on 
dynamic  grids,  PhD.  Dissertation,  Stanford 
University. 

Ramakrishnan  S,  Collis  SS.,  2002:  Variational  multi¬ 
scale  modeling  for  turbulence  control,  First  Flow 
Control  Conference,  AIAA  2002-3280. 

Roe  PL.,  1981:  Approximate  Riemann  solver,  parameters 
vectors  and  difference  schemes,  J.  Comp.  Phys.; 
43:357-371. 

Sagaut  P.,  2006:  Large  eddy  simulation  for 

incompressible  flows,  Third  Edition,  Springer- 
Verlag. 

Smagorinsky  J.,  1963:  General  circulation  experiments 
with  the  primitive  equations,  Mon.  Weather  Rev.; 
91:99-164. 

Tran  H,  Koobus  B,  and  Farhat  C.,  1998:  Numerical 
solution  of  vortex  dominated  flow  problems  with 
moving  grids,  AIAA  Paper  98-0766. 

Turkel  E.,  1987:  Preconditioned  methods  for  solving  the 
incompressible  and  low-speed  compressible 
equations,  J.  Comp.  Phys.  72:277-298. 

Van  Leer  B.,  1979:  Towards  the  ultimate  conservative 
difference  scheme  V:  a  second-order  sequel  to 
Godunov's  method,  J.  Comp.  Phys.;  32:361-370. 

Viozat  C.,  1997:  Implicit  upwind  schemes  for  low  Mach 
number  compressible  flows,  INRIA  Report  NO. 
3084. 


7 


